Two-qubit decoherence mechanisms revealed via quantum process tomography 
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We analyze the quantum process tomography (QPT) in the presence of decoherence, focusing on 
distinguishing local and non-local decoherence mechanisms for a two-partite system from experi- 
mental QPT data. In particular, we consider the \/iSWAP gate realized with superconducting phase 
qubits and calculate the QPT matrix \ in the presence of several local and non-local decoherence 
processes. We determine specific patterns of these decoherence processes, which can be used for a 
fast identification of the main decoherence mechanisms from an experimental x-matrix. 
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I. INTRODUCTION 

Quantum information processing is presently a focus 
of significant interest, since it shows promise to perform 
various computational and communication tasks which 
are difficult or impossible to perform by classical means 
[l| . A standard scheme of quantum information process- 
ing involves a sequence of unitary operations (gates) on 
single qubits or pairs of qubits. Due to coupling to en- 
vironment, the quantum-processor evolution suffers from 
decoherence, which introduces errors into quantum in- 
formation processing. Effects of decoherence on the de- 
sired quantum evolution can be characterized by a vari- 
ety of methods jointly called Quantum Process Tomogra- 
phy (QPT), such as the standard QPT [TJ, 2,9, ancilla- 
assisted process tomography (AAPT) [J, [f| m, nj , and 
direct characterization of quantum dynamics [8|, [9( . The 
standard QPT is simplest of the above methods in the 
sense that it can be performed with initial states being 
product states and local measurements of the final states. 

In recent years the QPT has been demonstrated exper- 
imentally in optics 0, M, [IIJI^JrI HBIli, NMR 

and for solid-state 

and three-qubit 
Experiments on the 



I, for ions in traps 
qubits [llJllJ^JSl^ne-qubii 
two-qubit [HlHlllMEIl 

HI HH systems have been studied 
QPT involving more than one qubit usually use the stan- 
dard QPT. The QPT experiments with the supercon- 
ducting phase qubits [H, [24|, [Hj , which are of the most 
interest for us here, have been also based on the standard 
QPT. In the present paper we also use the standard QPT. 

The QPT provides a very rich (complete) information 
on the performance of a quantum circuit. For N qubits 
the QPT matrix \ [H is generally characterized by 16^ 
real parameters; this number reduces to 16^ — 4^ param- 
eters if we limit ourselves by trace-preserving processes. 
Only 4 N — 1 of these parameters correspond to a unitary 
evolution, while the rest of them are due to decoherence. 
The problem of converting experimental QPT data into 
a characterization of decoherence processes is of signif- 
icant theoretical interest [H, S3, El HI H S3 • How- 
ever, there is still no good understanding of the relation 
between the x-matrix elements and decoherence param- 
eters, except for the case of one qubit. Recently, esti- 




mation of one-qubit decoherence parameters by a QPT 
method was discussed in Ref. [3(| for a specific decoher- 
ence model. In practice, decoherence models are often 
not known in advance, especially for systems contain- 
ing more than one qubit. An additional complication 
is that often two or more unknown mechanisms simul- 
taneously cause decoherence. In the present paper, we 
consider an approach to identify the decoherence mod- 
els from the form of the x-matrix provided by an ex- 
periment. For that we start with physically reasonable 
models of decoherence and analyze corresponding pat- 
terns in the x-matrix. If these patterns are sufficiently 
specific, then the main decoherence mechanisms can be 
identified from an experimental ^-matrix directly, with- 
out a complicated numerical analysis. As a particular 
example we consider the -s/iSWAP gate made of super- 
conducting phase qubits [H [H, Hj] and calculate the 
%-matrix in the presence of several local and non-local 
decoherence mechanisms, which can be anticipated for 
this system. We show that the patterns of significant 
elements of the ^-matrix are quite different for different 
decoherence mechanisms, that makes their identification 
relatively simple, even when two or more decoherence 
mechanisms simultaneously affect the system. 

The paper is organized as follows. In Sec. Ill Al we re- 
view the standard QPT for a generic system (with some 
formulas discussed in Appendix [A]) and then in Sec. Ill Bl 
we modify this formalism to make it more convenient 
for application to a bipartite system. Section IIIII is de- 
voted to a brief discussion of the Markovian decoherence 
and calculation of its contribution into the ^-matrix. In 
Sec. IIVI we introduce quantitative characteristics of the 
decoherence non-locality, which can be obtained from ex- 
perimental QPT data. Section|V]is the maj or part o f our 
paper, in which we analyze the two-qubit ViSWAP gate 
made of superconducting phase qubits. We start with the 
discussion in Sec. IV Al of an ideal ViSWAP gate, then in 
Sec. IV Bl we discuss several applicable models of local and 
non-local decoherence In Sec. IV Cl these models are used 
for the calculation of the ^-matrix of the trivial (iden- 
tity) two-qubit gate, then in Sec. IV Dl the %-matrix of 
the ViSWAP gate is calculated for the same decoherence 
models (it happens to have significant similarities with 
the identity gate case), and these results are discussed in 
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Sec. lVEl In Sec. IV Fl we analyze effects of decoherence on 
the x-matrix in the case of coupled but strongly detuned 
qubits. Section |Vl] is the brief conclusion. 

II. QPT BASICS 

A. QPT for a generic system 

According to quantum mechanics, a closed system un- 
dergoes a unitary evolution determined by the system 
Hamiltonian. However, usually quantum systems are 
coupled to environment, i.e., are open. The evolution of 
an open quantum system is described P, |35| by a com- 
pletely positive linear map C (a quantum operation): if 
initial density matrix of the system and environment at 
time t = is a product state, p° <g> p E , and full evolu- 
tion is described by Hamiltonian i?sE, then at time t the 
reduced density matrix of the system only is 

d-i 

P = £[P°], Pij = C HMPkU (1) 

kd=0 

where d is the dimension of the Hilbert space of the sys- 
tem, and the superoperator C has elements 

dm= E (ii'|c- iHBBt/B |fcA / >0* / |e-^ t / R |K / )*^. 

i' ,k' ,1' 

(2) 

with i,j,k,l denoting orthonormal basis states of the sys- 
tem and i' , k' , I' denoting the environment basis states. 

Besides the four-index quantity Cij^u h is convenient 
to introduce 0, [36| the d 2 x d 2 matrix C with the same 
components, but indexed in a different way: 

£{ij){kl) = £ij,kl, (3) 

where we use notation 

(ij)=di+j, (4) 

so that (ij) — 0, 1, . . . , d 2 — 1 (notice mnemonic rule that 
the d-nary representation of the number (ij) is "ij"). 
Now Eq. can be recast as 

P = Cp°, (5) 

in which p is a column vector obtained by placing the 
rows of p one after another and then transposing the 
result, p^) = p^. 

The standard QPT 0, || is based on a different but 
equivalent description of a quantum operation: 

d 2 -i 

p = C[p Q }= J2 XmnE m p°El (6) 
m,n— 

where E n are linearly independent operators (in d- 
dimensional Hilbert space) and x is a d 2 x d 2 Hermitian 



positive-semidefinite matrix, which fully characterizes 
the quantum operation. A quantum operation should 
not increase the trace of the density matrix, that leads 
to the condition [l|, Il2| 

d 2 -i 

J2 XmnElE m <I, (7) 
m,n— 

where I is the (d-dimensional) identity operator. (For 
operators the inequality A < B means that B — A is a 
positive operator.) For trace-preserving operations Eq. 
([7]) becomes an equality, while trace-decreasing opera- 
tions correspond to situations when the system leaves its 
Hilbert space or we consider a measurement with a par- 
ticular result. 

The QPT matrix x can be obtained from experimen- 
tal data in two steps, first calculating the matrix C and 
then converting it into the ^-matrix. To obtain C one 
needs to prepare d 2 linearly independent initial states p° 
(chosen out of experimental convenience), perform the 
evolution, and measure the resulting states p n using the 
quantum state tomography [l|, [3?|]- Using Eq. (0, we 
can write R = CRo, where R and i?o are d 2 xd 2 matrices 
constructed from p n and p® n as 

R(ij)n = (Pn)ij, {Ro)(ij)n = (Pn)y') ( 8 ) 

so that the nth column of R is p„, and similarly for Rq. 
Therefore, the matrix C can be obtained [2| as 

C = RRo\ (9) 

where the existence of Rq 1 is ensured by the linear inde- 
pendence of the states p° . 

Calculation of the x-matrix from C is the easiest when 
the operators E n used in the definition $Q form the "by- 
element" basis 

F m =\i){j\, (10) 

which we will call the "elementary basis" F n . This is 
because Eq. |T]) can be rewritten in a form similar to ^ , 

d 2 -i 

P= J ^nF m p°Fl (11) 

m,n— 

where d 2 x d 2 matrix J contains the same elements as £, 
but in a different order: 

J(ij}{M) = £(ik)(jl)- (12) 

Therefore, for the elementary basis, E n = F ni we obtain 
X = J, so that the ^-matrix consists of reordered ele- 
ments of C. Explicitly, this reordering is the following: 
1) each row of C is converted into a dx d matrix by se- 
quentially placing the strings of d elements below each 
other and 2) these matrices are placed from left to right, 
with a new row of matrices starting after each d steps. 
(Another reordering of the matrix elements of C is also 
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often used in the literature 0, H H, E3, HE El: the 



operator C, which is related to J as J, 



(ij)(kl) 



C 



(ji)(ik) 



[43j ■ Operators C and J are called Choi or Jamiolkowski 
operators. Both C and J are Hermitian and positive- 
definite.) 

To obtain x-matrix for a general operator basis E n , 
let us construct the d 2 x d 2 matrix E, whose nth col- 
umn contains all elements of the dxd matrix E n , so that 

,d 2 — 1 

F m E mn , and hence 
_1 exists because of 



Then E n = J2° m =o 

Fn = EL^".^" 1 )™- where E 
the linear independence of E n . In this way from Eq. (jllj) 



1 E„ 



we obtain 



x = E-\j(E-y. 



(13) 



use it in this paper. To perform the AAPT on a d-level 
system S, one needs a similar d-level ancillary system 
S' . The compound system is prepared in the maxi- 
mally entangled state 1$) = d -1 / 2 Ej=o K*)' t nen the 
quantum operation £ is applied to the system S only, 
and then the resulting density matrix of the compound 
system is measured by the quantum state tomography. 
It is easy to see that the resulting density matrix is 
(£ <g) 2)[|$)($|] = d- x J, where 1 is the identity map. 
In this way the matrix J is obtained directly and may 
later be converted into the %-matrix, as discussed above. 
In principle, other initial states can be also used for the 
AAPT, however the maximally entangled state |$) is the 
optimal one @. 



This expression simplifies in an important special case 
of mutually orthogonal operators E n , which satisfy equa- 
tion 



Tr(EtE n 



dS r , 



(14) 



where S nm is the Kronecker symbol and the convenient 
normalization factor d allows us to include the unity op- 
erator into the set E n (as well as products of Pauli matri- 
ces for multi-qubit systems) . Generalization to a different 
normalization is trivial - see below. In this special case 



E 



IT 1 * 

*,3=0 a (ij)f 



E, 



<*3> 



= Tv(ElE m ) = dd m 



(E'E) nm - 

i.e., E^E = dl (in other words, E j \fd is a unitary ma- 
trix), and therefore Eq. (fTB")) becomes 



X 



dr 2 E\jE. 



(15) 



In the case (|T4j) the calculation of the trace of the both 
sides of Eq. |(7J) results in the inequality 



Trx < 1, 



(16) 



which becomes the equality for a trace-preserving map. 

An important example of the orthogonal unitary- 
operator basis E n [satisfying Eq. (|14j) ] for a system of N 
qubits is the so-called Pauli basis, which consists of tensor 
products of N operators from the set {/, X, Y, Z}, where 
X, Y, Z are the Pauli operators. The modified Pauli basis 
with Y — > — iY is also used in the literature, for example 
in the QPT analysis for one and two qubits in Refs. [3,13]. 

Notice that if in Eq. (fT4")) the normalization factor d 
is replaced by an arbitrary number Q, then in Eq. (|15p 
the factor d~ 2 is replaced by Q~ 2 and Eq. (f!6|) becomes 
Tr% < d/Q. In particular, Q = 1 for the orthonormal 
basis F n introduced by Eq. (110|k in this case E = I, and 
therefore Eq. (fT5")) reduces to the previous result x = J. 

Several useful formulas for the ^-matrix are discussed 
in Appendix |A] Notice that the QPT calculation pro- 
cedure discussed above is slightly different and simpler 
than in Refs. P, @j; in particular, it involves an inversion 
of a d 2 x d 2 matrix [Eq. ([9])] instead of a pseudoinverse 
calculation for a d 4 x d 4 matrix. 

At the end of this subsection let us briefly discuss the 
idea of the AAPT @, H @, 0, 

even though we will not 



B. QPT for a bipartite system 

Now let us consider a bipartite system S consisting 
of subsystems Si and 52 with the Hilbert space dimen- 
sions di and c?2, respectively. Then the dimension of the 
Hilbert space of S is d — did2, and as the basis we can 
use the states 



\3i)\32) 



\3l32) 



(17) 



constructed out of orthonormal basis states in two sub- 
systems. We enumerate the states using slightly general- 
ized notation so that j — (j\j2) = d2.ji + j%. 

As discussed in the previous subsection, the %-matrix 
can be calculated by performing the quantum opera- 
tion on d 2 initial states p®. It is often convenient to 



use the product states, 



Prix 



(2) 



(where 



(711712) = d\ni + 1x2) with linearly independent sets of 
states for each subsystem. In this case the calculation of 
the matrix C via Eq. (JS|) may be simplified; however, this 
requires some modification [43j of Eq. ([9]). The reason is 
that the matrix Rq does not coincide with the Kronecker 



product .Rq 1 ' 1 <8> R [ Z> , as may be naively expected, but re- 
quires an additional permutation of rows. As the result, 
it is easier to calculate first the matrix C! = R[(R^ )~ 1 ® 
(i?Q 2 ^) -1 ], and then obtain C by permutation of columns, 
£m{hni 2 h) = An(iii 2 3*ii2>> where the four-number nota- 
tion in indices is the natural generalization of the no- 
tation (g]): {hjiiijz) = iidid 2 + j x d\ + i 2 d 2 + 32 and 
(i\ii3\3i) =i\d\d\ + i 2 did 2 + j\d 2 + j 2 . 



? (2) 



In particular, for a a two-qubit system with the initial 
states chosen as products of the states 0, 0, H| pn — 

pi? = |^n)<^n| With |^0> = |0), = |1>, 1^2) = 

(|0> + |1>)/V2, and |^3) = (|0) + we obtain 



(4 2) )- 



( 1 -(l+i)/2 (-l + »)/2 
-(l+»)/2 (— 1 + i)/2 1 
1 10 
V i -i 

(18) 

As discussed in the previous subsection, the calculation 
of the ^-matrix is the easiest when for the operator basis 
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E n we choose the elementary basis (fT0|) . Then \ = Jj 
where J is given by Eq. (|12p . However, for a bipartite 
system it is convenient to use the product of operator 
bases for each subsystem, 

E {nin2) =E$®E$, (19) 

and the product of the elementary bases for each subsys- 
tem is not the elementary basis ([TU)) because of different 
enumeration. Therefore, to simplify formulas for a bi- 
partite system, we have to somewhat modify the formu- 
las for the generic system. In particular, for the prod- 
uct Fn{ ® of the elementary bases (fTU|). we get 
X = J, where J(i 1 k 1 i 2 k 2 )(j 1 i l32 i 2 ) = _£(iii2jij 2 )(kik 2 hh}- 
(The relation between J and J is J(i 1 k 1 i 2 k 2 )(j 1 hj 2 i 2 ) — 

J(ili 2 k 1 k 2 )(j 1 j 2 l 1 l 2 ) ■) 

For the basis (|19[) which uses arbitrary orthogonal sub- 
system bases E„ and En^ , satisfying equations 

Tr(£«t£«) = di 5 nm) Ty(E^E^) = d 2 S nm , (20) 

the x-matrix can be expressed via J as 

X = d- 2 (E^ <g> E^)J{E^ (21) 

which is similar to Eq. (|15[) [a straightforward applica- 
tion of Eq. (TT5]) would not have the desired Kronecker- 
product form]. Notice that if subsystem bases satisfy 
the orthogonality condition (|2U)) . then the compound ba- 
sis (fT9|) satisfies the orthogonality condition ([14)) with 
the normalization factor d = did 2 - Therefore, Eq. (TIB|) 
remains valid, so that for a trace-preserving operation 

In particular, for a two-qubit system and the Pauli 
basis we have 



d = 4, E {1) = E {2) = 




III. MARKOVIAN DECOHERENCE 

A. General formalism 

An important special case of a general quantum evo- 
lution is the Markovian evolution, 

p = M[p], (23) 

where (M[p])ij = J2k~i=o ^ij,HPkl and the superoperator 
M is the generator of a quantum Markov semigroup [35| . 
The four-index representation of M can be converted into 
a d 2 x d 2 matrix M with JW/yw^ = MijM [similar to 
Eq. (3j)] , and then 

C = e Mt . (24) 



It is often convenient to separate the evolution gen- 
erator M = £ co h + L into the coherent part, £ C ohP = 
— (i/K) [H, p] with H being the Hamiltonian of the sys- 
tem, and the generator L of the incoherent evolution 
(decoherence) . In the matrix form we have 

M = L coh + L, (£ C oh)<y)(fci) = i(HijSik — Hi k Sji). (25) 

In the present paper we are interested in effects of deco- 
herence, and therefore we assume that the Hamiltonian 
H is known. Given the matrix £, which can be measured 
as discussed in Sec. HH the matrix M can in principle 
be extracted by solving Eq. (|24|) (t houg h the extraction 
procedure involves some subtleties) [22TI261. [2i| . Then L 
can be obtained from Eq. (f2"5)) . 

Following the QPT description ©, it is convenient to 
introduce a d 2 x d 2 matrix A defined by the equation 

d 2 -i 

L[P\= ^mnE m pEl (26) 

m,n=0 

with the same operator basis E n . The matrix A is Her- 
mitian, and for a trace-preserving map has 

Tr A = 0. (27) 

The matrix A is a counterpart of the ^-matrix and has 
many similar properties [43|. In particular, for a bipartite 
system and for a product basis E n satisfying Eqs. fTS"]) 
and (|20[) , the matrix A is given by an equation similar to 
Eq. (gU, 

A = d~ 2 (E^ <g> E^) v <g> E^), (28) 

where v^kx^k^ihhhh) = L ^mn){kik 2 i 1 i 2 )- Notice 
that A = v for the elementary product basis F^i <B> F^% ■ 

B. Weak decoherence 

The decoherence should be relatively weak for a prac- 
tical quantum information processing. In this case (for a 
sufficiently short time t) one can expand C up to the first 
order in L and obtain in the interaction representation 

£i*t-£i + r dre- L °°* T Le L °°* T , (29) 
Jo 

where C 1 = I is the c? 2 -dimensional identity matrix, and 
the interaction representation describes the evolution of 
p**(t) = e im ' n p{t) e~ lHt / h . 

Further simplification is possible for a very short time 
or when the secular approximation [44| is applicable [43| . 
Then the time-dependent factors in the integrand in Eq. 
(gg) can be omitted, yielding £ int = C 1 + Lt and 

X int = X 1 + At, (30) 

where x 1 is the process matrix for the identity map (see 
Appendix [X]). Unfortunately, the secular approximation 
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is usually applicable only when the subsystems (qubits 
or qudits) are uncoupled and there are no external fields, 
so that in the situations typical for quantum information 
processing (quantum gates) the simple equation (|30D is 
not applicable. Notice that the conversion between the 
Schrodinger and interaction representations for the x~ 
matrix is (see Appendix [A"| x = V\ lnt V\ where V is 
a unitary matrix with V nrn — Tr(E\e~ lHt / h E m )/d for 
the orthogonal basis E n satisfying Eq. (| 14[) . 



IV. 



CHARACTERISTICS OF NONLOCAL 
DECOHERENCE 



QPT provides a wealth of information: there are d 4 
independent real parameters in the matrix x ( or d 4 — d 2 
for a trace-preserving quantum operation) , and the num- 
ber of these parameters increases exponentially with the 
number of subsystems. However, the number of indepen- 
dent parameters for a multipartite system decreases dras- 
tically for local (independent) decoherence of the subsys- 
tems. In this section we discuss local decoherence of a 
bipartite system (generalization to a multipartite system 
|43| is rather straightforward). 



A process matrix \ is factorizable if and only if x — X- 
For x — X m a trace-preserving case when Trx' 1 ^ = 
Tr%( 2 ) = Tr% = 1, the matrices x^ 1 ' and x^ m Eq. (pTf]) 
necessarily coincide with x and X^ ■ 

If X 7^ Xi we can introduce a dimensionless parameter 
£nl characterizing non-locality of the decoherence: 



£nl = Tr|x - x|/Tr|x - Xidcail, 



(34) 



where Xidcai is the process matrix for the ideal coherent 
operation, which would occur in the absence of decoher- 
ence, and the absolute value of a matrix A is defined as 
| A | = V A^A, so that Tr\A\ is the so-called "trace norm" 
of A. Since decoherence yields a deviation of x from 
Xidcah the reasoning behind the definition (|3"4"|) is com- 
parison of matrices x ~ Xidcai and x ~ Xidcai, and char- 
acterization of their relative difference. For £nl -C 1 the 
factorization (|33[) is still a good approximation, while for 
£nl ~ 1 the decoherence is significantly non-local. Notice 
that the definition (|34p is meaningful only in the absence 
of Hamiltonian coupling between the subsystems. 



B. Coupled subsystems 



A. Uncoupled subsystems 

Let us start with assuming uncoupled subsystems, so 
that unitary evolution is local. If also decoherence is 
local, it is easy to show that for the product basis (| 19[) the 
X-matrix is the Kronecker product of the corresponding 
X-matrices for the subsystems, 



x unc = x (1) 



X 



(2) 



(31) 



In this case the number of independent parameters is 



d\ 



4 



(or n' = d\ + d\ - d{ 



J d\ in the trace- 



preserving case), which is much less than for a general 



X-matrix: n x 



d\d\ (or n' = d\d\ - d\d\). There is 



roughly a square-root decrease of complexity (ATth root 
decrease for an A^-partite system). In particular, for a 
two-qubit system n x = 32 and h' = 24 versus n x = 256 
and n' x = 240. 

In an experiment it is generally not known in ad- 
vance whether decoherence is local or not. Therefore, 
a quite important information can be obtained by check- 
ing whether or not a given x-matrix has the product form 
pip or, more generally, by quantifying the accuracy of 
the product-form approximation. 

Let us define the reduced x-matrices for subsystems as 



*«=Tr 2 x, X (2) =Tn 



X 



X 



(32) 



(in more detail, Xmlm = J2m 2 =oX(m im2 )(n im2 ) and sim- 
ilarly for x^)j and introduce 



x (1) ®x (2) - 



(33) 



In the case of Markovian evolution, the nonlocality of 
decoherence can be checked even in the presence of a 
coupling between the subsystems. We assume that the 
coupling is included into the (known) Hamiltonian H and 
that the generator of the incoherent evolution L [and 
hence the matrix A - see Eq. (f26|) ] can be extracted from 
experimental data. For the case of local decoherence the 
generators and of the subsystems decoherence 
contribute to L as [HI 

^ (iii2hh) (kik 2 hl 2 ) — ^ , \i > 1 j 1 )(k 1 l 1 )^^2k 2 ^]2h 

+ £ <*2j2)<fc2i2}^l fc l' 5 ?l'l' ( 35 ) 

and there is a simple relation 

A = A< 1 W (2) +X I{1) ®\ {2 \ (36) 

where x 1 ^ and x 1 ^ are the identity-map process ma- 
trices for the subsystems. 

Similar to the discussion above, we can introduce re- 
duced matrices -V 1 ) = Tr 2 A and \^ = TriA and their 
combination 



A = A«®x /(2) +X /(1) ®A^. 



(37) 



Also similarly, it can be shown [43| that a given matrix A 
has the local-decoherence form (j3"6")l if and only if A = A. 
In such a case A 1 - 1,2 ) = A' 1,2 -*, assuming trace-preserving 
operation with Tr A (1 ) = TrA< 2 ) = Tr A = 0. When A ^ 
A, the nonlocality of decoherence can be characterized by 
the dimensionless parameter 



4 L = Tr|A-A|/Tr|A|. 



(38) 
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Note that the nonlocality parameters e^ L and £nl are 

(1 2) 

invariant under a change of the bases E„ ' , which pre- 
serves orthogonality [Eq. (|20|) ]. Analysis J3| shows that 
£ nl ~ £ nl when decoherence is weak, the subsystems are 
uncoupled for coherent evolution, and either there is also 
no coherent evolution of the subsystems or the secular 
approximation holds. 

V. EFFECTS OF DECOHE RENCE 
MECHANISMS ON TWO-QUBIT ViSWAP GATE 

Even for only two qubits, the number of decoherence 
parameters in the x~ ma trix is quite big: in a trace- 
preserving case we have d 4 — 2d 2 + 1 = 225 parameters. 
This corresponds to the number of generally possible de- 
coherence processes. Obviously, interpretation of exper- 
imental ^-matrix data in such a case is quite difficult. 
However, instead of considering all general decoherence 
processes, it is meaningful to consider only physically 
reasonable mechanisms. Then by identifying specific fea- 
tures of these mechanisms in the ^-matrix and comparing 
with experimental data, it is possible to find the magni- 
tudes of various decoherence processes. 

In this section we consider the ViSWAP gate made of 
superconducting phase qubits [H, H3| and calculate the 
X-matrix assuming several plausible models of decoher- 
ence. We focus on identification of specific features of the 
X-matrix, which may serve as an evidence for a particular 
mechanism. In particular, we emphasize distinguishing 
local and non-local decoherence mechanisms. 

A. The ViSWAP gate 

The qubit states |0) and |1) of a superconducting phase 
qubit [32| are the ground and first excited states in the 
potential well. In this section (except Sec. IV F[) we as- 
sume that the two qubits are in exact resonance and use 
the rotating frame, which zeroes the Hamiltonians of the 
individual (uncoupled) qubits. Then the Hamiltonian of 
the capacitively coupled qubits in the rotating frame has 
the form [33|, |4f| 

H =(ftS/2)(|01><10| + |10><01|), (39) 

where S is the coupling strength (we assume that S is 
real). The Hamiltonian (|3"9"1) . which can be recast as 
H = (hS/4)(X <g> X + Y ® Y), is a special case of the 
exchange Hamiltonian, the so called XY Hamiltonian. 
It was extensively discussed in relation to quantum com- 
putation. Estimation of the exchange Hamiltonian by 
means of the QPT was discussed in Refs. [30l[3l|. 

The evolution of the two-qubit system is then described 
by the unitary operator 

U(t)= e - im ' n = |00)(00| + |11)(11| 

+ cos(Si/2)(|01)(01| + |10)(10|) 
-isin(Si/2)(|Ql)(10| + |10)(01|). (40) 



For a non-integer value of tS/2n the gate (|40| is an 
entangling gate and therefore, together with one-qubit 
gates, it is sufficient for quantum computatio n pal . In 
particular, U(w/S) provides the iSWAP gate [43], while 
U(it/2S) = ^vTswap is the ViSWAP gate For 
phase qubits the operation of the ViSWAP gate has been 
demonstrated experimentally [H, HH . 
We use the Pauli basis, 

E {nin2 ) =X ni ®X n2 , (41) 

where {X , X u X 2 , X 3 } = {I,X,Y,Z}, so that 
{E ,E 1} ...,E 15 } = {I ® 1,1 X,I ®Y,I ® Z,X ® 
Z (g) Z}. Note that 'nin^ is the base-4 represen- 
tation of (nxnz), e.g., Eg — X 2 ® Xi = Y ® X . The op- 
erators (|4Tj) satisfy the orthogonality condition (fl4|) with 
d = 4, so that Tr(E^E m ) = 45„ m . Any linear (Kraus) 
operator K in the 2-qubit Hilbert space can be repre- 
sented as 

15 

K = J2 k nEn (42) 

n=0 

with k n = Tr(E^K)/4. Correspondingly, any quantum 
operation of the form p = K p°K^ is described in the 
Pauli basis by the process matrix (see Appendix [A)) 

k m k* n . (43) 

In the Pauli basis Eqs. (J40J) — (J42J) (with K=U) yield 

U(t) = { [1 + cos(St/2)] sin(5t/2) (X <g> X 

+Y ®Y) + [l-cos{St/2)]Z®Z}/2, (44) 

and for t = ir/2S this becomes 

^v/Iswap^ I( 2 + V2)I®I-iV2(X® X + Y&Y) 
+(2-V2)Z®Z}/4. (45) 

The process matrices x f° r f ne gates (|44|) and (l45|) can 
be calculated using Eq. (|43[) . The process matrix Xideai 
for the perfect ViSWAP gate is shown in Fig. [1] (since x 
is Hcrmitian, the shown elements are symmetric about 
the main diagonal in the upper panel and antisymmetric 
in the lower panel). The nonzero elements of the matrix 
Xideai are 

Xoo - (3 + 2V2)/8, X15.15 = (3 - 2V2)/8, 

X55 = XlO.10 = X5.10 = Xl0,5 = X0.15 = Xl5,0 = 1/8, 

X05 = Xo,io = -X50 = -Xio,o = i(V2 + l)/8, 

Xl5,5 = X15.10 = -X5,15 = ~XlO,15 = «(V2 - l)/8. (46) 

An advantage of using the Pauli basis for the x-matrix 
of ViSWAP is that it results in a relatively small number 
of nonzero elements (8 real and 8 imaginary ones) out of 
the total number 256. For comparison, in the elementary 
basis ■[ 1*1*2} O'l J2 ] } the operator U v / igwAp has six terms 
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[see Eq. (|40|) with t = n/2S], resulting in 6 2 = 36 nonzero 
terms in the %-matrix. Notice that conversion of Xidcai 
from the Pauli basis to the modified Pauli basis (with 
Y — > —iY) would require sign change of six elements: 
Xio.o, Xio,5, Xio,i5, Xo,io, X5,io, and xi5,io (in general 
this conversion changes 138 out of 256 elements of the 
X-matrix: 18 elements change sign, 60 elements are mul- 
tiplied by i, and 60 elements are multiplied by —i). 




FIG. 1: The process matrix Xidcai for the perfect %/iSWAP 
gate in the Pauli basis. 



B. Models of decoherence 

In this subsection we consider several physically rea- 
sonable decoherence models for two phase qubits (all 
Markovian and trace-preserving) , including local deco- 
herence [H, [5(1 HH, EH [H, correlated dephasing 

f l56l and noisy coupling. In Refs. 0, [H^, HH, [52l 
154 I55L [El| these models have been mainly used to 
analyze two-qubit entanglement and Bell-inequality vio- 
lation, while in this paper we focus on their effect on the 
X-matrix of a quantum gate. Notice that estimation of 
one-qubit decoherence parameters by the QPT was dis- 
cussed in Ref. [301 ■ 



1. Local decoherence 

First, let us consider the model of local decoherence, 
described by the Bloch equations [H| for each qubit, so 
that in the rotating frame the density matrix of a sepa- 
rated qubit a (a = 1, 2) evolves as 



.(a) 

Pll = 
.(a) 

PlO = 



.(a) 

-Poo 



— _r( Q ) „( Q ) 1 -p(a) „(«) 

— 1 d Pll fl « 



'Poo 



-p\V/n 



(47) 



where and are the energy relaxation rates, so 



that T- 



(«) 



(lf } 



- 1 - ^ d t Tu)- 1 is the energy-relaxation 

time, and T 2 is the dephasing time (the Bloch equa- 
tions correspond to the secular approximation for a non- 
degenerate two-level system weakly coupled to a bath; 
r« Q ' ) = for a zero-temperature bath). 

Comparing Eq. (|47|) with the equation p a = L^ a 'p a , 
we obtain the one-qubit Markovian generators 



/ -r! tt) 













-1/T 2 (Q) 



















V rl Q) 








-r?° ; 



(48) 



while the two-qubit generator L\ oc of the local deco- 
herence is then given by Eq. (|35p . in which (ijkl) — 
8i + 4j + 2k + I (so that Hjkl' is the binary represen- 
tation of (ijkl)). 

Notice that the model of local decoherence involves 
two decoherence mechanisms: energy relaxation and pure 
dephasing. Correspondingly, X loc = Iq OCiE R + -^Ioc.pd- 
Technically, this splitting corresponds to representing de- 
phasing rates as sums of two terms, l/T^ — (rjj + 



T { u a) )/2 + T a , and then zeroing either T a or T d 



(a) 



In the presence of decoherence, the x-matrix typically 
acquires additional nonzero elements (in comparison with 
Xideai)- As shown below, the positions of the most signifi- 
cant extra elements of x may reveal the main mechanisms 
responsible for decoherence. 



2. Correlated dephasing 

Now let us consider two models of non-local decoher- 
ence, starting with the model of correlated pure dephas- 
ing. For a pair of coupled phase qubits, the correlated 
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dephasing can result from fluctuations of a common part 
of the magnetic field biasing qubits. We consider the sys- 
tem Hamiltonian H+HcT>(t), in which H is given by Eq. 
(f3"9"| . while the dephasing contribution is 



H CD {t)= a{*i(t)|io)<io| + <y a (t)|oi)<oi| 
+[*(*) + fc(t)]|ii><ii|}, 



(49) 



where Si (t) and <5 2 (t) are random but partially correlated 
frequency shifts for the two qubits. In the derivation of 
Eq. (|49|) we neglected noise-induced transitions between 
the levels, assuming that the noise intensity at the qubit 
frequency practically vanishes. 

Applying the standard method [H, [57], [H[ , we obtain 
the Markovian master equation for the average density 
matrix: 



P =-(i/n)[H,p] + LcD[p], 



^CD 



[p} = - 



r 2 p i 

r 2 pio 

TiP2o r_p 2 i 

T+P30 TiP3i 




(50) 



(51) 



where r± - V t + T 2 ± f , T a = f™{6 a (0)6 a (t)) dt (a = 
1,2), T = JZ° (5i(0)5 2 (t) + 5 2 (0)Si(t)) dt, and we have 
assumed (5 a (t)} = 0. The parameter of common dephas- 
ing r is zero in the case of uncorrelated (local) dephasing, 
while f = ±2v / TiT2 for full correlation/anticorrelation, 
so that the dimcnsionlcss correlation parameter is k = 
T /2y / Tir2, — 1 < k < 1. In the following subsections we 
will mainly focus on the case Ti = T 2 = Tpd- Notice that 
Eq. (|5ip is written in the computational basis | j) = | ji J2) 
with j = = 2j 1 + j 2 , so that j = 0, 1,2,3 corre- 

spond to ji]2 = 00, 01, 10, 11. In deriving Eqs. ([5TJ)l and 
(jtljl we have assumed r Q r c CD < 1 and SV C CD < 1, where 



is the correlation time of the frequency fluctuations. 
In discussion of the QPT it is very easy to get lost with 
different bases used in different equations. So we would 
like to repeat which bases do we use. In Eq. (|51|) [as 
well as in Eq. (|53[) below] we consider a two-qubit den- 
sity matrix, so this 4x4 matrix uses the two-qubit basis 
{|00), |01), 1 10) , |11)}. Then this equation is converted 
into the equation for the 16x16 matrix L, which uses the 
basis of 16 elements of the 2-qubit density matrix. The 
matrix C — e i L coh.+L)t ugeg same "by-element" basis 
as L. Finally, the matrix C is converted into the 16 x 16 
matrix \, for which we use the basis of product-Pauli op- 
erators. Somewhat differently, in the previous subsection 
Eq. (g7]) uses the one-qubit basis {|0), |1)} and Eq. Ijlgj) 
uses still one-qubit but 4-dimensional "by-element" ba- 
sis. Equation (|4"5|) is then converted into the equation for 
L, which uses the same 16-dimensional basis as above, 
and further procedures coincide. Notice that the bases 
discussed in this paragraph have nothing to do with the 
set of initial states discussed in Sec. [XT] [e.g., in the para- 
graph above Eq. (|18|) ]. which would be important in the 
experimental procedure. 



3. Noisy coupling 

The second non-local decoherence model we consider 
is the model of a noisy coupling. In the case of capaci- 
tively coupled phase qubits this model corresponds to a 
fluctuating coupling capacitance; a more practically im- 
portant case is when qubits are coupled via a tunable 
Josephson circuit, whose parameters may fluctuate. In 
the Hamiltonian (f3"9"|) we substitute S with S + s(t), as- 
suming (s(t)) = 0. Then following the same derivation as 
in the previous subsection, we obtain the master equation 



p=-{i/h)[H iP ]+L NC \p] ) 



(52) 



where 




-poi 

-2/311 + 2p 2 2 

2pi2 - 2p 2 i 

~P31 



~P02 

-2pi2 + 2p 2 l -Pl3 

2pil - 2p 22 -P23 
-P32 

(53) 

andT s = (1/4) f™ {s(0)s(t))dt. Inthederivat ion we have 
assumed T s t^ c < 1 and St^ c < 1, where t^ c is the 
correlation time of s(t). 



When the discussed above decoherence mechanisms ex- 
ist concurrently, the system state obviously evolves as 

p = -{i/K)[H,p] + (L loc + L C D+L NC )[p]- (54) 

By fitting experimental data with this model, it is pos- 
sible to find the corresponding best-fit decoherence rates 
quantitatively, and determine in this way if a particular 
decoherence mechanism is important or not. However, 
this is a rather laborious procedure. Another way to find 
out which decoherence mechanisms are important, is by 
checking characteristic features in the x-matrix, unique 
for a given mechanism. We will identify such features in 
the following subsections. 



C. Effects of decoherence on the identity gate 

Before studying the effects of decoherence on the x~ 
matrix of the ViSWAP gate (that will be done in the next 
subsection), let us consider decoherence for the identity 
gate, i.e., for the vanishing two-qubit Hamiltonian. Then 
C = e with the models for the decoherence generator 
L discussed above, and C can be converted into x m the 
way discussed in Sec. Ill Bl 

We are interested in effects of weak decoherence, cor- 
responding to sufficiently short gate-operation times. In 
this case the process matrix for the identity gate can be 
approximated [see Eqs. ((30|) and (|A8|) ] as 



X w X + At, 



(55) 



where A is determined by Eqs. (|28|) and (|22|) (we use 
the Pauli basis). The matrix A is a sum of contributions 
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from different decoherence mechanisms, which have the 
following explicit forms. 

For the local energy-relaxation mechanism, the 
nonzero matrix elements of A in the Pauli basis are 

Aoo = -2(rW + rV 2) ), 

An = A22 = r+\ A44 = Agg = r+\ 

Ao.3 = -^30 = T-\ Aq,12 = Al2,0 = T-\ 

A21 = -a 12 =ir (2) , A 84 = -A 48 = zr (1) , (56) 

where r£ a) = (I^ a) ±ri a) )/ 4 [notice a difference with the 
notation T± used in Eq. (|51[) ] . The contribution from the 
local pure-dephasing mechanism is a special case of the 
correlated dephasing which we discuss next. 

For the (correlated) pure dephasing the nonzero matrix 
elements of A are 

A o = -(ri + r 2 )/2, 

A33 = r 2 /2, Ai 2! i2 = 

^3,12 = ^12,3 — — -V),15 — — ^15,0 — r/4. (57) 

The absence of the correlation, f = 0, corresponds to 
the local pure dephasing; in this case the third line in 
Eq. l[5"7]) vanishes. 

For the noisy coupling the nonzero elements of A are 

Aoo = — r s , A55 = Aio,io = -^5,10 == Aio,5 

= -^0,15 = ^15,0 = r s /2. (58) 

All nonzero elements of x except for xooj are induced 
by decoherence. Because of the first-order approximation 
(|55p , the most significant additional elements of x are re_ 
lated approximately linearly to the nonzero elements of 
A (the second-order elements of x should be significantly 
smaller for a weak dephasing). Now, a very important 
observation is that the positions of the nonzero elements 
of A (excluding Aoo) in Eqs. (f56")) - (|5"5|) are different for 
different decoherence models, except for Ao^s and Ais^ 
appearing in both Eqs. (|57|) and (f58|) . Therefore, in the 
case of weak decoherence one can identify the consid- 
ered decoherence mechanisms simply by the positions of 
the most significant (first-order) elements of X- Another 
important observation is that effects of different decoher- 
ence parameters on elements of A are easily distinguish- 
able. In particular, for the local decoherence model [Eqs. 
([56]) and (f5"7"[) with T = 0], the decoherence in the first 
and second qubits is completely separated, affecting dif- 
ferent elements of A. Similarly, for each qubit the pure 
dephasing is separated from energy relaxation by affect- 
ing different elements of A, and the temperature for each 
qubit can be extracted from the ratio T^/T^, which is 
equal to the ratio of corresponding elements of A in Eq. 
([56")) . The correlation factor in the pure dephasing model 
can be extracted via the relative height of elements in 
the second and third lines of Eq. |[5T|) (two positive ele- 
ments in the third line should be used, since the negative 
elements are also involved in the noisy coupling model). 





ERt>o 


ERt— 


LPD 


LDt>o 


LDt=o 


CD 


CCD 


NC 


A 


13 


13 


3 


15 


15 


7 


7 


7 


V 


32 


23 


12 


32 


23 


12 


10 


16 


X 


64 


64 


4 


64 


64 


8 


8 


8 


J 


36 


25 


16 


36 


25 


16 


16 


20 



TABLE I: The number of nonzero elements in the matrices A, 
v, X) an d J for several decoherence models: energy relaxation 
(ER) with arbitrary or zero temperature T, local pure dephas- 
ing (LPD), local decoherence (LD) which includes both pure 
dephasing and energy relaxation (with T > or T = 0) , corre- 
lated pure dephasing (CD) with < |k| < 1, completely cor- 
related/anticorrelated pure dephasing (CCD) with k = ±1, 
and noisy coupling (NC). We use the Pauli basis for the ma- 
trices A and Xj while in the elementary basis they would be 
equal to the matrices v and J, correspondingly. The identity 
gate is assumed for matrices x an d J. 



The clear separation of effects allows us to estimate the 
relative values of the decoherence parameters simply by 
the relative values of the corresponding elements in the 
X-matrix. 

Such a simple analysis is possible, to a large extent, 
because we use the Pauli basis. The matrix A in the 
Pauli basis has a relatively small number of non-zero ele- 
ments. The first row in Table[T]shows this number for the 
considered decoherence models. For comparison, in the 
elementary operator basis { |*i«2) (J1J2 1 } we have A = v 
[see Eq. (|28[)]. and then the significantly larger number 
of non-zero elements is given by the second row in Table 
[I] [the number of nonzero elements of v equals that of L] . 
Since Eq. (|55[) is only an approximation, the number of 
non-zero elements of the matrix x is larger than that for 
A; it is shown in the third row of Table Q] for the Pauli 
basis. For the elementary basis [then x = J, see Eq. (|21[) ] 
this number is shown in the fourth row and is typically 
significantly larger (except for the model of energy relax- 
ation). This illustrates convenience of the Pauli basis. 

The nonlocality parameters £nl and e^ L introduced 
by Eqs. (|34[) and l[55|) approximately coincide in the 
weak-decoherence case (|55[) ; e^ L is time- independent, 
while £nl slowly changes with time. Calculation of e^L 
gives the following results for the considered decoherence 
models. For the local decoherence involving both en- 
ergy relaxation and pure dephasing, we obtain e^ L = 0, 
as should be expected. For the model of correlated 
pure dephasing with Ti = T2, the nonlocality param- 
eter e^ L = 2 1 k; |/(1 + \/i + k 2 ) depends on the corre- 
lation factor k, so that e^ L « |k| for |k| 1 and 
e^ L = 2\/2 - 2k 0.83 for k = ±1. Finally, for a noisy 
coupling e^ L = 2-^2— 1 w 1.83. As expected, e^ L is of 
the order of 1 for a strongly nonlocal decoherence. 
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D. Effects of decoherence on the ViSWAP gate 

Now let us co nsider t he effects of decoherence on the x 
matrix of the \ASWAP gate. We calculate \ via Eq. (12T 



from the evolution equation C = eS Lcoh+L " ,7T ^ 2S , where 
i co h is given by Eqs. (f25) and (f59"| . and 16 x 16 matrix 
L depends on the decoherence model fSec. IV Bj) . 

In the important case of weak decoherence the first- 
order approximation ([29| leads to the linear relation be- 
tween the decoherence contribution x — Xideai and the 
decoherence generator L (the evolution time is fixed). 
Therefore, since the decoherence generators for various 
mechanisms simply add up [see Eq. (|54|) ] , their contribu- 
tions into the \ matrix are approximately additive for a 
weak decoherence, 



X ~ Xideai + Sxioc + SXCD + #XNC, 



(59) 



and we can consider them separately. 

Figures [2H1] discussed below show the numerical results 
for the x-matrix of the v iS WAP gate in the presence of 
the decoherence mechanisms considered in Sec. IV Bl A 
comparison of Fig. Q]with Figs. EMU shows that the effect 
of decoherence on the x-matrix j s t modify the values of 
the nonzero elements of Xideai and generally to add extra 
nonzero elements. Below we identify patterns of extra 
elements specific for each considered decoherence mecha- 
nism, which allow for a fast, though tentative, identifica- 
tion of these mechanisms. For this purpose it is sufficient 
to consider only the largest extra elements (most signifi- 
cant out of the first-order in decoherence elements). An 
important observation (see below) is that for the con- 
sidered decoherence models the positions of the largest 
extra elements of x coincide with the positions of ele- 
ments of A discussed in the previous subsection. This 
makes the analysis of the x-matrix for the ViSWAP gate 
rather similar to the analysis in the absence of unitary 
evolution (Sec lVCj) . except for the noisy-coupling model 
for which there are no extra elements in x- Consider now 
the effects of decoherence models in more detail. 



1. Local decoherence 

We focus on parameter values typical for experiments 
with the superconducting phase qubits, and therefore as- 
sume zero-temperature case (which is a ve ry good ap - 
proximation for the experiments [23l. l24l. l25l l32l. [33l. |34| ) . 
We also assume the same local-decoherence parameters 



for both qubits, so that T 



(i) 



.(2) 



(1) 



T, 



(2) 



T 2 , and rl 1} 



F (2) 



0. 



Figure [2] shows the process matrix of the ViSWAP gate 
in the presence of local decoherence. For this example we 
have chosen the coupling S/2-k — 20 MHz (which is in be- 
tween the coupling values of experiments [33[ and (34[) 
and the decoherence parameters 7^=90 ns and T 2 — 60 ns, 
which are also more or less typical for the superconduct- 




FIG. 2: The process matrix % (i n the Pauli basis) for the 
ViSWAP gate in the presence of local decoherence for S/2tt = 
20 MHz, Ti = 90 ns, and T 2 = 60 ns. The matrix elements 
marked by the long arrows, X33 — Xi2,i2, are the features of 
the pure dephasing, while the elements marked by the short 
arrows (xn = X22 = X44 = Xss and X03 = X30 = Xo,i2 = Xi2,o 
in Rex and X21 = -X12 = X84 = ~X48 in Imx) are the 
features of the energy relaxation. 



ing phase qubits (much longer relaxation times have been 
achieved recently [H, H3|). 

The local decoherence includes two mechanisms: en- 
ergy relaxation (with the rate 1/T\) and pure dephasing 
(with the rate r PD = 1/T 2 - l/2T x ). As follows from 
the results presented below, the relative strength of these 
two mechanisms can be easily estimated by inspection of 
the extra elements of x (compared to Xideai) in Fig. [2] 
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The elements marked by the long arrows are due to pure 
dephasing, while the elements marked by the short ar- 
rows are due to the energy relaxation. By comparing the 
hight of the elements of x marked by the long and short 
arrows, one can crudely estimate the relative strength of 
these two mechanisms. 

The largest extra elements for the energy-relaxation 
model (short arrows in Fig. ^ in the first order in 1/Ti 
are the following: 



Xll = X22 = X44 = Xi 
X03 = X0.12 = X30 = Xl2,0 
X21 = X84 = -X12 = -X48 = i- 



n + 2V2 _ 0.37 

16ST1 ~ ~sT\ ' 

?r(2 + V2) _ 0.34 
32ST X ~ ST\' 
.7T + 2V2 _ .0.37 
WSTt ~ l ~ST~i 



(60) 



They are at the same positions as the elements of A, 
Eq. ([56]) (the remaining element Aoo is at the location 
of the main V iS WAP peak, while we consider only extra 
elements of x~ matrix). 

We emphasize that the ^-matrix also contains many 
other first-order in 1/Ti elements (in contrast to the 
unity-gate case considered in the previous subsection); 
however, they happen to have relatively small absolute 
values. The largest of them are imaginary: X46 — X4,n — 
X13.6 = Xi3.il = Mr/(16V2STi) « O.Ui/Sn; this gives 
8 elements together with the corresponding Hermitian 
conjugated elements. There are also 4 elements of mag- 
nitude 0.06/STi and 8 elements with the absolute value 
0.02/STi. 

The rest of the extra elements of x are of a higher order 
in 1/jSTi, and therefore much smaller than the first-order 
elements for ST\ 3> 1. For instance, X33 — X1242 = 
X3,i2 — Xi2,3 ~ (ti' 2 /64)(S'Ti) _2 (we show these elements 
explicitly because they are located at the same positions 
as the most significant elements for the model of pure 
dephasing discussed below). 

For the model of local pure dephasing the largest (not 
all) first-order in Tpo elements are 

(3^ + 2)r PD 0.71 r PD 

X33 = X12.12 ~ 7777; ~ a , (61) 



16S 



s 



and they are again at the same positions as the elements 
of A in Eq. (|57p. These elements are shown by the long ar- 
rows in Fig. [2 Since for Fig. [5] we assumed Tpd = 1/Ti, 
the height of these elements is comparable to (approxi- 
mately twice larger than) the height of the main extra 
elements due to the energy relaxation. The other (much 
smaller) first-order in Tpd elements are discussed below, 
combined with the more general case of correlated pure 
dephasing, which we consider next. 



X3,12 — Xl2,3 



r/2rpD- Now x generally contains eight extra elements, 
all of them real. These eight elements are also present in 
the first order in Tpd/S; however, only four of them are 
relatively large: 

[3tt + 2 + (tt - 2)k] r PD 

X33 — Xl2,12 ~ ^Qg 

w (0.71 + 0.07k) r PD /5, (62a) 
[tt - 2 + (3tt + 2)n] r PD 
165 

« (0.07 + 0.71k) r PD /5, (62b) 

while the other four elements are much smaller, X66 — 
X99 = ~X69 = ~X96 ~ (tt - 2)(1 - K)r PD /(165) « 
0.07(1 - K)r PD /5. Notice that the larger elements ([62)) 
are again at the positions of the elements of A in Eq. 
([ST)) , and these positions are all different from those for 
the energy relaxation. For a weak correlation, k<1, the 
elements X3,i2 and Xi2,3 [see Eq. (|62b[) ] become small, re- 
covering the result (|61j) for the local dephasing. 



Figure [3] shows the x-matrix of the V iS WAP gate af- 
fected by the partially correlated dephasing, k = 0.5. 
The two diagonal extra elements, X33 and Xi2,i2, are at 
the positions shown by the long arrows in Fig. (2[ and 
their values are almost independent of k. In contrast, 
the off-diagonal elements X3,i2 and Xi2,3, marked by the 
arrows in Fig. (3[ strongly depend on the correlation k, 
so that their magnitudes are comparable to the values of 
X33 and Xi2,i2 on ly f° r a significant correlation, |«| ~ 1 
[see Eq. (j62b[) ]. This clearly suggests the way to check 
decoherence due to fluctuating common magnetic field in 
an experiment with phase qubits. 

The x-matrix has especially simple form in the case of 
the completely correlated dephasing, k — 1, so that T_ = 
0, and r + = 4FpD in Eq. ([5l[) . Then the exact solution 
gives the following nonzero elements of x- The elements 
which are at the positions of the nonzero elements of 
Xidcal become 



t u 



X:„ 



1 1 
1 1 



-ig- 
-ig- 



(63) 



V 11 ig- ig- f- ) 

where m, n = 0, 5, 10, 15, f± = 2 + ± 2V2j d , g± = 
\/27d±l, and 7^ = e ~ 7rrpE> / 2S . The extra nonzero matrix 
elements are 



X33 — X3,12 — Xl2,3 — Xl2,12 — (1 



-27rr PD /S 



)/8, (64) 



so that all of them are equal [as for the first-order result 
([6"2"| with k = 1], Notice that for a partially correlated 
dephasing with < K < 1 the exact solution gives X33 = 
Xi2,i2 > X342 = Xi2,3 > (as in Fig.©. 



2. Correlated dephasing 



3. Noisy coupling 



Let us consider the effects of correlated pure dephasing 
for Ti — T2 = TpD and arbitrary correlation factor k = 



In contrast to the previous models, the matrix A for 
the noisy-coupling decoherence [Eq. ([B"8"[) ] has non-zero 
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0.51 





.... j^m^mffimm 






"■pgr 






FIG. 3: The x-matrix (in the Pauli basis) for the ViSWAP 
gate in the presence of correlated pure dephasing for S/2ir = 
20 MHz, Tpd = (90 ns) _1 , and k = 0.5. Significant matrix 
elements shown by the arrows, X3,i2 = Xi2,3, indicate signifi- 
cant correlation of the two-qubit dephasing, \k\ ~ 1. 



FIG. 4: The x-matrix (in the Pauli basis) for the ViSWAP 
gate in the presence of a noisy coupling for S/2n = 20 MHz 
and r s = (90 ns) _ . The specific feature of the noisy cou- 
pling is the increase of the elements shown by the arrows in 
comparison with their ideal value (v2— l)/8. 



elements only at the positions, for which the matrix Xidcai 
is also non-zero. As a result (not quite trivial), the noisy 
coupling does not produce extra elements in the x-matrix 
[neither in the first-order approximation (|29[) nor in the 
exact solution]. The exact solution gives the following 
non-zero elements of x : 





f 3 


f 2V2 7c 


ih+ 


ih + 


1 \ 


1 




— ih+ 


1 


1 


—ih_ 


8 






1 


1 






V 


1 


ih- 


ih- 3 


- 2a/2 7c / 



where m,n = 0,5,10,15, h± — v27c ± 7c , and j c = 
e -7rr 3 /25^ The x _ matrix for S / 27T = 20 MHz and T s = 

1/(90 ns) is shown in Fig. 2J 

Despite the noisy coupling does not produce ex- 
tra elements in the x-matrix, this model also has its 
unique feature. Let us consider the imaginary elements 
Xs,is, Xw,is, Xib,5, and xi5,io, shown by the arrows in 
Fig. |U which all have the same absolute value h-/8. 
From the above formula it is easy to see that this value 
is larger than the ideal value (y/2 — l)/8 ~ 0.052 (unless 
decoherence is very strong, T s /S > 0.77), with the max- 
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imum of 0.094 at T s /S — 0.22. In all other considered 
models the absolute value of these matrix elements de- 
creases in comparison with the ideal case [see Figs. HKS] 
and Eq. (|()3"|) ]. so their increase is a unique feature of the 
noisy coupling. 

Notice that the absence of this evidence does not ex- 
clude the possibility of noisy coupling, since the increase 
of the elements marked by the arrows in Fig. 3] may be 
compensated by the their decrease due to other deco- 
herence mechanisms. Generally, the fast identification of 
decoherence models by their unique features should serve 
only as a preliminary step, while the accurate quantifica- 
tion of the decoherence mechanisms requires a numerical 
best-fit procedure. 

E. Discussion 

As observed and discussed above, for the considered 
decoherence models the positions of the largest extra ele- 
ments in the x-matrix of the V iS WAP gate coincide with 
the positions of nonzero elements of A. Moreover, a com- 
parison of Eqs. ([HD and with Eqs. and (j5T)) 
shows that for the largest extra elements Xmn ~ A m „T ff , 
where r g = n/2S is the gate operation time [this state- 
ment is not correct for Eq. (|62b[) in the case |k| <C 1, but 
then the elements are small anyway]. This fact is not 
trivial and deserves a discussion. 

Let us consider an arbitrary two-qubit entangling gate 
described by a Hamiltonian H with a characteristic fre- 
quency S. In the presence of a weak Markovian deco- 
herence L the first-order approximation (|2Q[) gives the 
evolution map 

C(r g ) = e L ^ T " + f " e L "^ {T ^ T) L e L "^ T dr, (66) 
Jo 

where r g is the gate operation time. This matrix can be 
transformed into Xi giving the corresponding separation 
of terms \ = Xidcai + For a very short time so that 
Sr g <C 1, the exponential factors in Eq. (|6"6")1 are close to 
one, and therefore 

6x~\r g , (67) 

seemingly explaining the fact observed above. The prob- 
lem, however, is that a strongly entangling gate (such 
as ViSWAP) necessarily operates in the different regime, 
Sr g > 1, for which the approximation (|67|) is not valid. 

We have numerically checked the relative error of the 
approximation (|o7|) for the %/ iS WAP gate by calculat- 
ing the dimensionless parameter e = Tr\5x — Ar 9 |/Tr|(5x|, 
introduced similar to Eq. As expected, we have ob- 

tained e ~ 1 (e.g., e = 0.54 for the energy relaxation), 
confirming that the approximation (|67[) is invalid. How- 
ever, the inaccuracy e happens to be mainly due to a large 
number of small non-zero elements in Sx, while for the 
largest extra elements (where Xideai is zero) the approx- 
imation (|67p unexpectedly survives. The origin of this 



fact is still a puzzle for us. Nevertheless, the same use- 
ful property may hopefully be also valid for some other 
quantum gates and decoherence models. 

F. Strongly detuned qubits 

In Sec. IV CI we have calculated the x-matrix for the 
identity gate, which is realized when the two qubits are 
uncoupled. However, in many experimental realizations 
(e.g., for capacitively coupled phase qubits [33|, [3J|) the 
qubits arc permanently coupled, and effective uncoupling 
is produced by strong detuning of the qubits: |Aw| ^> \S\, 
where Aw = uJ q b,i — Uqbfl is the detuning. In such ex- 
periment it is natural to extract decoherence parameters 
from the x-matrices in both situations: resonant qubits 
(ViSWAP) and strongly detuned qubits (that gives a 
more straightforward access to decoherence parameters) . 
In this subsection we analyze the x-matrix in the case of 
strongly detuned qubits. 

To deal with detuned qubits we need to introduce a 
rotating frame, which rotates with different frequencies 
for different qubits. For this frame it is preferable to 
use the actual eigenfrequencies (shifted due to the level 
repulsion), even though we still use the basis of uncou- 
pled states ( 1 01) , 1 10} , etc.). This produces Hamiltonian 
[which replaces Eq. (|3"9f ] 

H = (h/2)(Au- Aw)(|l)(l| ®I- I® |1)(1|) 
+ (ftS/2)(e- lAlit |01) (10| + e l ^ M \\Q) (Oil), (68) 

where Aw = -J (Aw) 2 + S 2 x sgn( Aw) is the difference of 
eigenfrequencies (we define it with the same sign as for 
Aw), and the first term is due to the level repulsion. We 
assume strong detuning, which means |Aw| ~ |Aw| S> 
| S |; then the first term in (|68)l is small since Aw — Aw w 
5 2 /2Aw. 

For strongly detuned qubits, Eq. ([55]) for the identity 
gate x-matrix is somewhat modified, and at sufficiently 
small time t the x-matrix can be approximated as 

X « X 1 + $X C + At, (69) 

where x 1 is for the ideal identity gate, the small cor- 
rection 5x c comes from the Hamiltonian (|68p . and the 
decoherence A-matrix is also somewhat modified (as dis- 
cussed below). 

The correction Sx c oscillates in time with frequency 
Aw. After averaging over these fast oscillations we get 
(in the first order in 5*) four non-zero terms: 

SXo9 = *Xm = 'Woe = -SXso = *5/4Aw. (70) 

If we do not average over time, then these terms should 
be multiplied by (1 — cosAwt); also there will be four 
more terms with zero average: <5xo5 = ^Xo 10 = ~^xio = 
— Sxio o = *<S'sin(Awt)/4Aw. Notice that 8x° is small 
only in the rotating frame based on the eigenfrequencies. 
If, for example, the rotating frame uses the unperturbed 
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qubit frequencies, then both qubits will be slowly rotat- 
ing about the z-axis of the Bloch sphere, that will eventu- 
ally produce large terms Sxo 3 , 6x30, ^Xo,i2) S X c i2,0' S Xi2,3> 
6x3 i2- In an experiment the choice of the rotating frame 
corresponds to the choice of the reference microwave fre- 
quencies. 

Now let us discuss contributions of the dephasing pro- 
cesses to the matrix A. It can be shown that the contribu- 
tion from the energy relaxation in qubits is still given by 
Eq. ([56]) (as for the identity gate); however, the up/down 

rates r„ and for the two qubits (a — 1,2) can now 
be different from these rates for the qubits in resonance. 
The difference is because these rates are proportional 
to the Fourier components of the bath spectral density 
at the qubit frequencies ±u> q i, , Q , and therefore changing 
qubit frequencies may noticeably change the rates. For 
pure dephasing the contribution to the A-matrix is still 
given by Eq. ([57]) (as for the identity gate) without any 
changes. 

In contrast, the contribution to the A-matrix due to 
the noisy coupling S + s(t) significantly differs from 
Eq. ([55]) for the identity gate. This comes from a 
significant change of Eq. (|53[) when we take into ac- 
count the new Hamiltonian (|68p. First, the rate T s = 
(1/4) (s(0)s(t))dt should now be replaced by T' s = 
(1/4) (s(0)s(t)) cos(Au)t) dt. This would lead to a 
negligible change if the correlation time °f the noise 
s(i) is short: r c NC < l/|Aw|; however T' s < T s for a long 
correlation time: t^ c 3> l/|AcD|. The second change 
in Eq. ([53]) is that we should delete the term P21 in the 
second row and similarly the term p\2 in the third row 
(this simple change happens in the secular approxima- 
tion, F^ <C |A(D|). As a result of the changes, the con- 
tribution to the A-matrix due to noisy coupling has now 
the following nonzero elements: 

^00 = — r' s , Ao,i5 = ^15,0 = r s /2, 

A55 = Aio,io = A5,io = Aio,5 = ^66 = A99 = r s /4, 
A 69 = A 96 = -T'JA. (71) 

Notice four new positions of non-zero elements in this 
matrix compared to the identity gate case (|58[) . 

It is important to mention that positions of non-zero 
elements of Sx c and non-zero elements of A for vari- 
ous decoherence models are still all different (except el- 
ement Aoo, which is non-zero even in x 1 > an d elements 
Ao,i5)Ai5,o which appear in both correlated dephasing 
and noisy coupling). Therefore, measuring x-matrix ex- 
perimentally for strongly detuned qubits gives an easy 
way to find main decoherence mechanisms and quantify 
their parameters. Notice that measuring x-matrix ([S5]) 
for several times t gives a more accurate value for A by 
a simple least-square method and also allows for check- 
ing the linearity in time, which essentially checks that 
decoherence is Markovian. 



VI. CONCLUSION 

In this paper we have discussed the effects of deco- 
herence on the quantum process tomography of a quan- 
tum gate. In particular, we have introduced (Sec. IIV[) 
dimensionless parameters, obtainable from experimental 
QPT results, which characterize nonlocality of decoher- 
ence. As an important practical example (Sec. W\ . we 
have analyzed the process matrix x f° r the two-qubit 
%/ iS WAP gate in the presence of several local and non- 
local decoherence mechanisms, typical for superconduct- 
ing phase qubits. Besides presenting explicit results for 
the x-matrix in the presence of decoherence (using the 
Pauli basis), we have focused on finding specific patterns 
for each decoherence model. These patterns may be used 
for a fast identification of the most important decoher- 
ence mechanisms in an experiment, that is an alternative 
to the laborious procedure of numerical best-fitting of ex- 
perimental x-matrix. Somewhat unexpectedly, we have 
found that these patterns for the considered decoherence 
models are to a large extent the same for the ViSWAP 
and identity gates. In future it is interesting to study 
whether or not our fast-identification approach can be 
applied to other quantum gates and decoherence mecha- 
nisms. 
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APPENDIX A: SOME PROPERTIES OF THE 
PROCESS MATRIX x 

1. Let us consider the change of the matrix x under a 
linear transformation of the basis {£",1} — > {E' n }. From 

d 2 — 1 

Eq. ©, using the substitution E n = J2 n ,n'=o^n'nE n ,, 
where V is the d 2 x d 2 transformation matrix (while E n 

are dx d matrices), we obtain p = £m,nLo x' mn E' m P° E n 
with 

x = v x vK (Al) 

If the both bases are orthogonal, so that Tr(E n E m ) = 
QS nm and Ti(E$E'J = Q'S nm , then V nm = 
Tr(E'^ E m ) I Q' '; in this case \/Q'/QV is a unitary ma- 
trix. 

2. Let us consider the transformation of the x-matrix 
under unitary transformations of the initial and final 
states, p — » p' — UpW and p° — > p' = Uop°U^, where 
U and Uq are unitary operators. From Eq. © we obtain 

P' = Y^Ly=oXmnE' m p' E$, where E' n = UE n U ] , so the 
extra evolution of states corresponds to the transforma- 
tion of the basis E n . If the operators E n are orthogonal, 
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Tr(£t£ m ) = Q Snm> then also Tr(E$E'J = Q5 nm and, 
as follows from the previous paragraph, 

d 2 -l 

P'= E xLnEmp'oEl (A2) 

m,n— 

with given by Eq. (lAljl . in which V is now a uni- 
tary matrix with the elements V nm = Tr(E^E' m )/Q = 
Ti{EiUE m ul)IQ. 

3. Let us obtain the ^-matrix for an evolution p = 
Kp a K if with an arbitrary linear operator K . The most 
important special case is the unitary evolution (then K 
is unitary): however, in general K is an arbitrary Kraus 
operator [jij . Representing K in the operator basis E n as 

d 2 -i 

K=J2 k n E n (A3) 
)i=0 

and comparing the evolution equation with Eq. ([5]), we 
obtain 

Xmn = k m k n . (A4) 

Notice that for the orthogonal basis, Tr(E^E m ) = QS nm , 
the coefficients in Eq. (j4"!2|) can be calculated as 

k n = Tr{ElK)/Q. (A5) 



4. As a simple example, let us consider the process 
matrix x 1 f° r the identity map. In this case Ju^/jn — 
^■\ij){ki) = ^ikSji, and from Eq. (fT3|) we obtain 

xlin = X] (- E! ~ 1 )m(«>(- E_1 )n<ji)- ( A6 ) 

For the orthogonal basis, T^E^-E^) = Q5 nm , from Eqs. 
([A!)) and (|A"5j) with K = I,we find 

Xm„ = 0" 2 (TrS m )*TrS„. (A7) 

This equation further simplifies when Tr_E„ = for all n 
except for, say, n — (as in the case of the Pauli basis). 
Then the basis orthogonality yields Eq = y^Q/dl, and 
Eq. (|A7|) becomes 

Xmn = (d/Q)S m0 S n0 . (A8) 
For the usual normalization Q = d, it becomes xL = 
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